Effect of electron-lattice interaction on the phase separation in strongly correlated 
electron systems with two types of charge carriers 
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The effect of electron-lattice interaction is studied for a strongly correlated electron system de- 
scribed by the two-band Hubbard model. A two-fold effect of electron-lattice interaction is taken 
into account: in non-diagonal terms, it changes the effective bandwidth, whereas in diagonal terms, 
it shifts the positions of the bands and the chemical potential. It is shown that this interaction 
significantly affects the doping range corresponding to the electronic phase separation and can even 
lead to a jump-like transition between states with different values of strains. 
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I. INTRODUCTION 

A typical feature of the strongly correlated electron 
systems is the formation of inhomogeneous states^ 2 -. The 
nature of such inhomogeneities based on the electron cor- 
relations. However, their specific manifestations can in- 
clude the effects of different degrees of freedom existing 
in the solids: spin, charge, orbital, and lattice^. A strong 
electron-lattice coupling plays a fundamental role in such 
actively studied systems as high-temperature supercon- 
ductors, manganites, cobaltites, and other related mate- 
rials^. An important characteristics of all these mate- 
rials is a complicated electronic structure involving two 
or more conduction bands. The simplest model allowing 
appropriate description of electron correlation effects is 
Hubbard model and its multiband generalization-^^. 

In Ref. 13 it was demonstrated that existence of two 
bands in the Hubbard model gives rise to the existence 
of the electronic phase separation even in the absence 
of any type of additional factors (spin, charge, orbital, 
etc.). However, these factors are necessary to give a 
realistic description of the particular physical systems. 
For example, taking into account spin and orbital vari- 
ables allows a detailed picture of the phase diagram for 
manganites--' 11 . Including in the model the possibility 
of the spin-state transitions gives an explanation of the 
phase separation in cobaltite s 12 ' 13 . An important and 
vast field of research is a problem of existence of inho- 
mogeneities in high-temperature superconductors, espe- 
cially in cuprates 5 . In this field, the multiband Hub- 
bard model provides some insight in the properties of 
the cuprate superconductors. 

Electron-lattice coupling was incorporated in the 
multiband Hubbard model to describe electronic struc- 
ture of cuprates and manganite o 14 ' 15 : 16 ] 17 . In these pa- 
pers, the main emphasis was put on the polaron effects or 
the influence of the electron-phonon interaction on elec- 
tron pairing. Here we use a similar type of the electron- 
lattice coupling to analyze its effect on the formation 
of the inhomogeneous states within two-band Hubbard 
model. 



II. MODEL AND A QUALITATIVE ANALYSIS 

Let us consider a strongly correlated electron system 
with two types of charge carriers, a and b, interact- 
ing with lattice strains u (static dispersionless phonons). 
Here, we limit ourselves to the case of small strains, 
\u\ -C 1, when the theory of elasticity is applicable. 

The Hamiltonian of such a system can be written as 



H — H e i + Hjj + H e i_ph + H p h- 



(1) 



Here, H e i corresponds to the energy of charge carries 
without taking into account the interaction between them 
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where a. and a- la (b\ and b- la ) are the creation and an- 
nihilation operators for a (b) electrons at site i with spin 
projection <r, (...) means the summation over the nearest 
neighbors, t a and t b are the corresponding hopping inte- 
grals, AE is the energy shift between a and b bands, and 
n? CT = a^dio- and n\ a = b\ a b\ a are the number operators 
for a and b electrons, respectively. 

The Hjj term corresponding to the on-site Coulomb 
repulsion of charge carriers has the form 



ni s + U b nln\ s )+^ 



Uab 

2 



(3) 



where U a , U a , and U a b are the energies of Coulomb re- 
pulsion between two a, two 6, and one a and one b 
electrons at one lattice site, respectively, and a means 
the spin projection with the sign opposite to that of a. 
We assume that the on-site Coulomb repulsion is large, 
U a ,Ub,U a b 3> zt (z is the number of nearest-neighbor 
sites, z = 6 for the simple cubic lattice considered in this 
paper) and U a ~ Ub ~ U a b ~ U. 

The electron-lattice interaction can be chosen in the 
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following form 

Hei- P h= ( Aa L a j CT "j + X ' b h b i* u i) ( 4 ) 

+ ^ (Aalii?!?, + AfcUiTl^) , (5) 
i.a 

where A, A' and A 0i f> are corresponding constants of the 
electron-phonon interaction and u\ are the local distor- 
tions corresponding to site i. Thus, we take into account 
the effect of lattice strains both on the on-site electron 
energy and intersite charge transfer. 

We approximate the phonon self-energy term H p h as 
an elastic energy of the system depending on distortions 
at different sites 

^"oiast = y^u?, (6) 

i 

where K is the elastic modulus. 

We analyze the problem in adiabatic approximation 
considering phonons as a classical static elastic field. To 
find a self-consistent solution to the problem, we first 
perform averaging of the Hamiltonian over the electronic 
degrees of freedom. From the condition of the energy 
minimum with respect to strains, d(H (ui))/du\ — 0, we 
obtain 

u i- K i {<) 

where i and j are the nearest-neighbor sites. 

Using Eq. (JT]), we can present the effective electron 
Hamiltonian as 

H * = - E - A ^ a l a >° + (* b - A ' Si ) b l b ^ 

- [AS + (A Q - A 6 )u,] n{ a + H v 

icr 

- 5> - * a ui) « + nt) + f E^i' ( 8 ) 

i,cr i 

here \x is the chemical potential. 

Hamiltonian Eq. (JSj) clearly demonstrates that the ef- 
fect of electron-lattice interaction is two-fold. This inter- 
action in non-diagonal terms changes the effective band- 
width whereas in diagonal terms, it shifts the positions 
of the bands and the chemical potential. In the earlier 
analysi s 9 ' 18 , we have shown that the qualitative features 
of the phase diagram for the two-band Hubbard model 
are mainly determined by two dimcnsionless parameters: 
the ratio of the bandwidths and the relative positions of 
the bands. Thus, to construct a minimal model capturing 
the main physical effects of electron-lattice interaction, it 
is sufficient to keep only A and A&. In addition, we put 
AE = to emphasize the effect of band shift related only 
to the electron-lattice interaction. As a result, we get 

_ _ Az(a^a jfT ) + \ b {n%) 



and 

H °x = - E [(*° - Afi ') a l a i° + thh l h -A + H v 

+ A b ui< CT -5> k + ni) + f E • ( 10 ) 

icx i.a i 

Starting from Hamiltonian (|10[) . we can point out the 
main qualitative effects of the electron-lattice interaction 
in the two-band model. In the absence of doping, n — 
0, bands a and b are empty and their centers coincide. 
With the growth of n the wider band a begins to be filled 
up from the bottom. The band filling is accompanied 
by strain Ui — —\z(a\ a a- ia )/K. The average (a\ a a- ja ) 
is proportional to the hopping probability and thus is 
positive. The strain Ui is positive if A < and negative if 
A > 0. From Eq. (flQ|) . it is easy to see that at any sign of 
A the bandwidth increases with the strain. At a certain 
doping level, the chemical potential attains the bottom 
of the narrower band b and this second band starts to be 
filled up. In this range of doping, we have two types of 
the electrons and the energy of the system depends on n 
in a rather complicated manner due to electron-electron 
correlations. As we have shown earlier—, such situation 
is favorable for the phase separation even in the absence 
of the electron-lattice interaction. 

However, the characteristic feature of the system un- 
der study is the dependence of the effective band shift 
AE c g — — XbUi on the strain and, hence, on the doping, 
according to Eq. ([9]). The sign of the shift depends on 
the signs of A and Af,. A simple analysis of Eqs. © and 
(fTr]|) shows that for the same signs of A and Af, the value 
of AE c ff is negative and the sign of u\ remains the same 
at any n. If the signs of A and Af, are different, the strain 
can change its sign at some doping level. The change of 
the sign of the strain results in the change of the sign of 
the effective band shift AE c q. 

The dependence of AE e R on doping can give rise to 
a more sophisticated behavior of the system. If at some 
doping level n* , the narrower band crosses the bottom 
of the wider band, that is \%n* / K > zt a , then it could 
be favorable to have almost all electrons in the b band. 
So, there appears a competition between two states with 
different values of the strain. It suggests the possibility 
of a transition between these two states, which can have 
a jump-like form. 

In the next section, we present a quantitative analysis 
of the possible situations. 

III. MEAN FIELD APPROACH 

In the limit of strong electron correlations, U — > oo, we 
can describe the evolution of the band structure with the 
change of n following the method presented in Ref. 
We introduce one-particle Green functions for a and b 
electrons. For band a, we have 

G aa (i- j, t) = -i{fa ia (t)al(0)}, (11) 
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where T is the time-ordering operator. The similar ex- 
pression can be written for band b. The equation of mo- 
tion for the one-particle Green function with Hamiltonian 
(flT)]) includes two-particle Green functions of the form 

GavM'ii ~ J> t) = -i(Taia(t)n b icr ,(t)al(0)) . (12) 

In the considered limit of strong on-site Coulomb repul- 
sion, the presence of two electrons at the same site is un- 
favorable, and the two-particle Green function is of the 
order of l/U. The equation of motion for two-particle 
Green functions includes the three-particle terms com- 
ing from the commutator of a\ a {t) or bi a {t) with the U 
terms of Hamiltonian (fTOf . which are of the order of l/U 2 
and so on. In these equations, following the Hubbard I 
approach^, we neglect the terms of the order of l/U 2 
and use the following decoupling in the Green functions 

(Ta i+mff (*)nU*)aj CT (0)) - (nl,)(fa i+m „(t)al(0)). As 
a result, we derive a closed system for the one- and two- 
particle Green functions 9 '- 9 . This system can be solved 
in a conventional manner by passing from the time-space 
(t, r) to the frequency-momentum (oj, k) representation. 
We limit ourselves to consideration of the case when the 
total number of electrons per site does not exceed unity, 
n = n a + n b < 1. The upper Hubbard sub-bands are 
empty, and we can proceed to the limit U — > oo. In this 
case, the one-particle Green functions are independent of 
U and can be written in the frequency-momentum rep- 
resentation as^ 

G a<T (k,w) = — - 



9a 



G ba (k,u>) = 



9aaW a (l 
9ba 



\u/t a )((k) 



+ (1 - X b u - g bcr WbC{k) ' 



(13) 



where we put ui — u assuming a homogeneous strain and 
introduce the following notation w a = zt a , a — a or 6, 



9n 



1 



(14) 



n" = (fifo) is the average number of electrons per 
site in state (a,cr), and £(k) is the spectral func- 
tion depending on the lattice symmetry. In the 
considered case of the simple cubic lattice, £(k) = 
— [cos(fc 1 <i) + cos(k 2 d) + cos(fc 3 d)] /3, d is the lattice 
constant. In the main approximation in l/U, the mag- 
netic ordering does not appear and we can assume that 



1 /2. Below we omit spin indices. 



IV. RESULTS AND DISCUSSION 
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Equations (fl3|) and (fl4j) demonstrate that the filling 
of the bands depends on the strain u and the number of 
electrons in one band depends on that in another band. 
Using the expression for the density of states p a (E) = 
— 7r _1 Im J G a (k, E+i0)d 3 k/ (2tt) 3 , we get expressions for 
the numbers of electrons in bands a and b 



z 2g a n 



fi 



g a w a (l-\u/t a ) 



, n b = 2g b n 






. gbWb 



(15) 




FIG. 1: (Color online) Evolution of the occupation numbers 
of the a and b bands with doping n: n a (red) solid line and n b 
(blue) dashed line. The parameters are X/w a ~ 0.8, \b/w a = 
2.0, K — lGuia and u>b/w a = 0.25. 



where 



n (fi') = J dE'po(E') 



(16) 



and p (E') = J d 3 kS(E' - C(k))/(27r) 3 is the density of 

states for free electrons (with the energy normalized by 
unity, \E\ < 1). The chemical potential /i in Eq. (fT5|) can 
be found from the equality 



n = n a (p) + n b (p,). 



(17) 



Equations (fT5|) include average strain u, which itself de- 
pends on n b and an average (al a aj a ). The latter can be 
expressed in terms of the Green function as 



-iG afT (i - j, 
du>d 3 ]i 



(2tt) 



-o) 

G aCT (k, uj)e 



(18) 



ik(i— j)+za;0 



Thus, we have a system of four equations ©, (fT5|) . and 
(fT7)l for finding u, n a , n b , and p. These equations are 
solved together with the relationships (fTB")) . p^|) . (|16p . 
and (fT8|) . 

The dependence of n a and n b on doping n is illustrated 
in Fig. [1] At low doping, the bottom of b band lies far 
above the bottom of the wider band a, and there exist 
only a electrons. The filling of the bands gives rise to a 
non-zero strain u, see Eq. @, and hence to the band shift 
AE e s . The plots u(n) and AE e g(n) are shown in Fig. [5] 
When A and Ab have the same signs the band b shifts 
downward (AE c g > 0). Thus, with the increase of n, 
the chemical potential crosses the bottom of the b band 
and b electrons appear in the system. Due to electron- 
electron correlations, the effective width of a band starts 



4 



-0.04 



IS 




).0 0.2 



0.6 n 0.8 




FIG. 2: (Color online) The dependence of u (solid (red) curve) 
and A.Eefi = —XbU (dashed (blue) curve) on n. The parame- 
ters are the same as in FigQ] 
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FIG. 3: (Color online) The energy of the system vs doping 
level n. Solid (red) curve corresponds to the homogeneous 
state, whereas the dot-dash (green) curve is the energy of the 
phase-separated state. The parameters are the same as in 
Figd] 



to decrease. This band narrowing and increase of AE c g 
leads to decreasing of the number of a electrons, and at 
some doping level, the charge carriers of type a disappear 
in the system, see Fig. [1] 

The energy of the system in the homogeneous state, 

Shorn / p a {E')E'dE\ 

is the sum of electron energies in all filled bands. After 
straightforward calculations, we can write £hom in the 
form 



Shorn = 2ff„iy a I 1 - 7^ ) £o 



+ 2g 2 w b e 



to 



gbWb 



9aW a (1 - Au/ta) 

Ku 2 
+ \bun b H — , 



(19) 



FIG. 4: (Color online) Lattice strain as a function of doping. 
The parameters are \/w a = —1-2, \b/w a = 2, K = lQw a and 
Wb/vda = 0.25. Jump-like transition between two states with 
different values of lattice distortions occurs at n — n* . Solid 
(red) lines correspond to the states with minimum energy, 
whereas dashed (blue) lines correspond to metastable states. 
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FIG. 5: (Color online) The energy of the system vs doping 
level n. Solid (red) curve corresponds to the homogeneous 
state, whereas the dot-dash (green) curve is the energy of the 
phase-separated state. The dashed (blue) curves correspond 
to the energies of metastable states. The parameters are the 
same as in Fig [3] The kink in the curve -Ehom(^) corresponds 
to the jump-like change in the lattice distortion. 



where 



e (ti') = J dE'E'p (E'). 



(20) 



The dependence of Eh om (n) is shown in Fig. [3]by solid 
line. We see that within a certain n range the system can 
have a negative compressibility, d 2 E\ wm /dn 2 < 0, which 
means a possibility for the charge carriers to form two 
phases with different electron concentrations^. 

The values of parameters chosen to plot Figs. [T][3] cor- 
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respond to a continuous evolution of the strain with dop- 
ing. However, as it was mentioned in Section II, one 
could also expect a jump-like transition between states 
with different values of the strains at certain values of 
parameters. Such a situation is illustrated in Figs. [4] and 
[5j In Fig. IU we can see that at some n — n* the strain u 
exhibits a stepwise transition with the change of the sign. 
In the vicinity of n* , there exist two competing states, 
A and B, corresponding to two solutions of the system 
of equations ©, (fl"5|) . and (fTTJ) for u, n a , n b , and fi. In 
the state A, we have n b <C n a or = 0, whereas in the 
B state, b electrons are prevailing. The energies of these 
states coincide at n = n* and at n > n* , state A has the 
a higher energy than state B. The minimum energy of 
the homogeneous state is shown in Fig. [5] by solid line 
and the energies of the metastable states are depicted by 
dashed lines. The change in the type of state corresponds 
to the kink in -Ehom(?i) curve at n = n*. 

The existence of two competing states in some range 
of parameters can be illustrated in the following way. 
Let us study Hamiltonian (|10p where the strain it = u is 
considered as an independent parameter. We analyze this 
Hamiltonian in the way similar to that described above. 
Namely, at each given u, we solve the system of equations 
(fl~5|) and (fTTf for n a , n b , and //, and then find the system 
energy per lattice site Ehomiu). The optimum value u 
is then determined by minimization of Eh om (u). The 
numerical analysis shows that the function E% om (u) has 
one or two minima depending on model parameters. The 
functions Ef,, om (u) calculated for two sets of parameters 
A, Xb, K and Wb at different doping levels n are shown 
in Figs [6] At u — u corresponding to the minimum of 
Ehomiu), we have 



dE h0 



On 



1 

N 



dH, 



du 



1 = ^ z ( a L a j°") + ^ br>b + = , 

(21) 

and we come back to Eq. © for u. 

In the general case, it is difficult to find explicit condi- 
tions for the existence of the jump-like transition. Here, 
we analyze the important particular case of A = 0. Let 
us consider the function E'{u) = dEh om /du, which now 
has a form E'(u) = Ku + Xbn b (u). The E'(u) curves 
calculated at different values of doping are illustrated in 
Fig. El For large negative strain, the b band lies far be- 



low the a band, and we have n a = and n 



Thus, 



E'(u) = Ku + Xbn linearly grows with u up to the point 
u = iti, when the chemical potential fi reaches the bot- 
tom of the effective a band —g a wa = — (1 — n)w a , and 
a electrons appear in the system. Using second equation 
in the system Eq. (| 1 5[) with nb — n, and [i = —(1 — n)w a , 
we obtain the following expression for u\\ 



-1 



Xi, 



(71\ ( IX 



, (22) 



where — 1 < fio(n) < 1 is the function inverse to n (u'), 
Eq. (flT)]) . This function is shown in the inset in Fig. [7] 
For large positive it, when band b lies above band a, we 
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FIG. 6: The energies of the homogeneous state depending on 
the strain u calculated for \/w a — 0.8 (upper panel), and for 
X/w a = —1.2 (lower panel) at different values of doping n. 
Other parameters are the same for both panels; Xb/iv a = 2, 
K = 16w a , and Wb/w a = 0.25. In the case of X/w a = —1-2, 
the function Ehom(u) has two minima within a certain doping 
range. 



have n a = n and n b = 0, and consequently E'{u) = Ku 
linearly decreases with it till the point it2 > u±, where b 
electrons appear. Acting in a similar way, we find 



1 

U 2 = — 



(1 - n)w b + (l - w a u a , 
It is clearly seen from Fig. [71 that if 



(23) 



E'(ui)=Ku 1 + X b n>0, E'(u 2 ) = Ku 2 < 0, (24) 

then the function E'(u) has three zeros, and the energy 
E(u), in turn, has two minima. Using Eqs. ((^)) and (|23[) . 
and taking the equality signs in relations (|24]) , we get the 
estimate for the minimum value of At,, at which the jump- 
like transition can occur 
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FIG. 7: The dependence of E'(u) on the strain u at different 
values of doping n. The parameters are A = 0, Xb/w a = 6, 
K = 16w a , and wt/w a = 0.25. The values of ui(n) and 
U2(n) > ui(n) are marked by circles. In the inset, the function 
(M){n), inverse to the function (|16[) is shown. 



where n<i is found from the equation 

n 2 \ 1 - n 2 Wb 



Mo 



2-712 



1 - n 2 /2 w a 



(26) 



The value of Aj* decreases with wj,. For very small ratio 
Wb/wa <C 1, we have /io [712/(2 — "-2)] ~ 0, that is, n 2 ~ 
2/3, and 



E'{u) can continue to grow (decrease) above (below) U\ 
(u 2 ), and, consequently, A£ found from these conditions 
overestimates its value. 

For A — the jump-like transition occurs for rela- 
tively large values of A&. For example, at K — ld>w a 
and Wb — 0.25w a , from Eq. (|25p one obtains A£ ~ 3.6w a - 
The numerical analysis shows, however, that even small 
negative (if Ab > 0) A sufficiently reduces the threshold 
value of A^. For example, at A = — 0.4w a , the jump-like 
behavior arises starting from A£ w 2.0w a (K = 16w a and 
Wb = 0.25w a ). Thus, different signs of A and At, favor the 
existence of such transition. 



V. CONCLUSIONS 

The electron-lattice interaction plays an important role 
in the systems with strongly correlated electrons affecting 
their behavior with doping. The interaction of electrons 
with the lattice distortions results, first, in the hopping 
probability and, second, in the relative shifts of the elec- 
tronic bands. We analyzed the problem in the framework 
of the two-band Hubbard model. We demonstrated that 
if the electron-lattice interaction is strong enough, there 
appears a competition between states with different val- 
ues of strains and the transition between these states can 
occur in a jump-like manner. We also showed that the 
electron-lattice interaction produces a pronounced effect 
on the conditions of the electronic phase separation since 
it influences the value of the bandwidth ratio and the 
relative positions of the bands. 



w a K 



(27) 



Note that the function E(u) can still have two minima 
when conditions (1241) are not met, since the derivative 
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